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We isolate the two-step mechanism involving a real intermediate photon from the one-step mech- 
anism involving a virtual photon for the trident process in a constant crossed field. The two-step 
process is shown to agree with an integration over polarised sub-processes. At low to moderate 
quantum non-linearity parameter, the one-step process is found to be suppressed. When the pa- 
rameter is large, the two decay channels are comparable if the field dimensions are not much greater 
than the formation length. 
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I. INTRODUCTION 

Partly due to experiments that have measured them, 
and partly due to theoretical proposals to observe them, 
higher-order quantum-electrodynamical processes in ex- 
ternal fields have recently gained much attention in the 
literature. Theoretical results for two-photon non-linear 
Compton scattering in a pulsed laser field [l], have 
shown in particle spectra a much richer physics of higher- 
order processes compared to tree-level versions. Recent 
attention has also been focused on the trident process 
in an external field, which is essentially lowest-order 
fcrmion-seeded pair creation, — > + e + e~, where 
e + represents a positron and e~ an electron. Part of the 
trident process was measured in the landmark E-144 ex- 
periment at SLAC US], 

which still more than a decade 
later is being analysed by theorists [5|, |6| despite higher- 
order processes, including trident, having first been stud- 
ied some years before 0-Q (a review of strong-field ef- 
fects in quantum electrodynamics (QED) can be found 
in iSp). 

In light of several plans to construct the next- 
generation of high- intensity lasers [l2j|, there has been 
much activity in attempting to simulate relativistic plas- 
mas that include strong- field QED effects l_3j. Due to 
their complexity and the current lack of a consistent 
framework for including classical and quantum effects 
alongside one another, approximations must be employed 
(although some non-perturbative strong-field QED sim- 
ulation methods for systems with fewer particles are also 
being currently developed [13, EBI)- The current letter is 
motivated on the one hand by the need to justify approx- 
imating higher-order QED processes by chains of tree- 
level processes in simulation-based approaches, and on 
the other by an enquiry into the physics of the trident 
process in an external field. This study complements the 
numerical approach of @ that analysed the weakly non- 
linear regime in E-144 using a monochromatic plane wave 
background modified to take into account finite interac- 
tion time, the lucid general theoretical outline of @ and 



the expression for the total rate in a constant crossed 
field neglecting exchange terms derived in @, [1] ■ By de- 
riving an analytical expression for the trident process in 
a constant crossed field, we will separate off in an unam- 
biguous way, the two-step process, measured in E-144 in 
a laser pulse, of a real photon produced via non-linear 
Compton scattering decaying into an electron-positron 
pair (e^ — > +7, 7 — > e + e~, where 7 represents a pho- 
ton). Moreover, we will show that the two-step process is 
exactly given by a sum over intermediate photon polar- 
isation of each tree-level sub-process integrated over the 
photon lightfront momentum and that for small quantum 
non-linearity parameter, the one-step process involving a 
virtual photon suppresses the total rate. Furthermore, we 
compare the relative importance of the two decay chan- 
nels and comment on the measurability. 

The paper is organised as follows. We begin by high- 
lighting important points in the derivation of the tri- 
dent process in a constant crossed field, relegating tech- 
nical albeit standard steps to the appendix. The two- 
step contribution is analysed and compared to combin- 
ing tree-level rates and then the remaining, nominatively 
"one-step" contribution is analysed and compared to the 
Weizsacker- Williams approximation. The total creation 
probability is then studied, the measurability of each pro- 
cess commented on, the results discussed and the paper 
concluded. 



II. PROBABILITY DERIVATION OUTLINE 

A diagram of the considered trident process is given 
in Fig. [1] where double lines indicate fermions dressed 
in the external field, which has a vector potential A^(tp), 
phase tp — xx and wavevector x, satisfying xA = x 2 = 
0. Following standard Feynman rules (see e.g. [Ta|). 
in a system of units c = h = 1 with the fine-structure 
constant a = e 2 , for positron charge and mass e > 0, m, 
the scattering matrix for this trident process is given by: 
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where the electron in, electron out and positron out wave- 
functions in the field of a plane wave ip, ijj, ip + are given 
by Volkov states [I?}, G lllJ (x — y) is the photon propa- 
gator and (p2 Pz) refers to an exchange of P2 and p^ 

in the first term of Sfi and Sfi = Sfi — Sfi. The sec- 
ond term must be subtracted due to exchange symme- 
try as the two outgoing electrons are indistinguishable 
(Pauli's principle). To avoid ambiguity arising from the 
exchange term, differential rates will be in p', p_ and p + , 
which refer to the scattered electron, created electron and 
positron momentum respectively. 




FIG. 1. The Feynman diagram for one term of the trident 
process in a plane wave (the other is given by the substitution 

P2 P3)- 



Let us fix the co-ordinate system by defining x^ = 
0,0,1)^, = A(¥>)(0,1,(L0)'\ Then focusing 

on just Sfi (the calculation for Sfi is analogous), us- 
ing the definition of the objects in Eq. ([I]) and Fourier- 
transforming both vertices x and y, one arrives at fl8j |: 



S~f\ = (2TT) 2 aJdrdsS {i ' ) (U)r l "{r) 



is 



k=k 



A M 00,(2) 



where II = p2 + pj, + Pi — p\ — (r + s)?c, k' is the 
photon wavevector, k^ = p\ — p2 + rx and T' i (r) and 
A M (s) are functions of variables at the first and second 
vertices respectively. It has been shown that the Fourier- 
transformation variables r and s are equivalent to the 
number of external-field pho tons, when the background 
is an infinite plane wave |10|. 

A constant crossed field background A** (ip) — a^ip 
is interesting, first because many integrals can be per- 
formed analytically facilitating physical interpretation, 
second, that integration is computationally sufficiently 
cheap that rates could feasibly be added to simula- 
tions and third that predictions in a constant crossed 
field are often a good approximation to in an arbitrary 
background field. When one considers that a general 
strong-field QED process can depend on four gauge- and 
relativistic- invariants 1 1 Oil 



T 



e 2 p^p„ 
m 2 (xp) 2 ' 
e 2 F fW F^ 
4m 4 : 



ey/\p„F^\ 



e 2 F*F^ 

Q 

y ~ 4m 4 ; 



(3) 



where T^" and F^ v are the energy-momentum and 
Faraday tensors and £ and x the classical and quantum 
non-linearity parameters, it is a common argument [Tol | 
that if £ 3> 1 (equivalent to process formation lengths 
being much smaller than the external field wavelength) 



the external field can be considered constant during 
the process, and if F, Q -C x 2 ,l, then probabilities P 
are well-approximated by those in a constant crossed 
field P{x,F,G) ~ P{Xi^i 0)- The classically non-linear 
regime £ 3> 1 is fulfilled by the most intense lasers [20| . 
as are F,G <^1. 

The probability of the trident process can be calcu- 
lated by performing the trace average over spin states 
(achieved using the package Feyncalc (2l| ) and integrat- 
ing over the outgoing degrees of freedom (a factor 1/2 
removes double-counting from identical final particles), 
P = (l/4)nJ=2[V / d 3 Pj/(2 7 r) 3 ]tr|S' /i | 2 , where V is the 
system volume. When the momentum-conserving delta- 
function in Eq. is squared, a factor in the denomina- 
tor of a formation phase length, A.(p + , is generated: 



6{r + s - (r' + s')) 



r-\-s—r' -\-s' 



A^+/2tt, (4) 



where A</? + = J dip+, (p± = tp x ±(p y is assumed finite and 
tp z = zx, z G {x,y}. The formation phase length can be 
related to particle momenta by calculating the position of 
the (real) saddle-point in the phase of Sfi, tp* + = <p x + p> y 
and then associating A<^!j_ = J dtp* + analogous to tree- 
level calculations [lOj, explained in more detail in Sec. 
IA 61 Replacing phase lengths Aip with coherence inter- 
vals Aip* is a key part of this approach and they will be 
used interchangeably. 

As the rate is proportional to \Sfi\ 2 = \Sfi + S~Ji\ 2 , we 
note interference between exchange terms arises. In Sec. 
IA 81 it is argued that this interference is negligible when 
the field dimensions are much larger than the formation 
length. This is the only part neglected as we proceed 
with P ?a (~P* + P)/2. Moreover, p2 «-» P3 is a symmetry 
of the remaining integrand, permitting us to define P = 
I 5 . 



Parts of tr are independent of outgoing momenta 

in the 1- (electric-field) direction, to differing degrees. 
Using the relationship between ip* + and p\ 3 , an integral 
J appears in P with the form: 



J 



1 

2^2 



dtp i dp_ 



dr 



- r F(r) 



(r + ie) 



(5) 



where F(r) £ C°°. A crucial step is how to deal with the 
integration over the photon propagator. As commented 
in jfjj, using the Sokhotsky-Plemelj formula (22J: 



dr 



F(r) 
r ±is 



TinF(0) + P / dr 



F(r) 



(6) 



where V refers to taking the Cauchy principal value of the 
integral, P might thought to be split into real and virtual 
parts, for which the photon is on-shell (k 2 = 0) and off- 
shell, corresponding to the first and second terms in Eq. 
© respectively. However, for a constant crossed field at 
least, these two terms correspond only to a 5 function 



3 



and a principal value part. Using Eq. ([6]), performing 
the principal values first, J can be shown to be equal to 



(7) 



= 2|F(0)| 2 / dip* + dcp*_ e(-<p*_) (8) 



& = « J*p* + jT* + ^(-0-2^(0) 

(9) 
(10) 



J, 



(i) 



d - - / 



d<p*l_ / dr 



|F(r) -F(0)| 2 



where = 2KeJ^\ Re is the real part and #(•) is 

the Heaviside theta function. J^' and J^ 1 ' 1 denote the 
direct one-step and one-step-two-step interference inte- 
grals. Since F most often has a maximum at r = 
for the most important dynamical regions, it can also be 
noted that the interference term is effectively negative. 
We return to this point in the discussion of the one- 
step process. Recognising that if ip* x ,Lp* y € ¥>] then 

2- 1 / <kp* + J dip*_6{-y*_) is equivalent to dip* f£ dip* y 

and forms a two-step process. Although there are 
ostensibly two phase regions tp* and ip*, in a constant 
crossed field, the two-step phase factor is trivially related 
to the total phase region through 



, , [*' , , (A^) 2 

fa Jfo 



(11) 



where Aip = ip — ipo is the total phase difference be- 
tween the beginning and end of the process in the ex- 
ternal field, thereby measuring its extent. What re- 
mains in P are terms proportional to 2 _1 j dip*^, equal 
to f^dip' = Aip, which we deem accordingly a one-step 
process. The Heaviside theta function in the two-step 
process preserves causality, ensuring that pair-creation 
from a photon occurs after photon emission from non- 
linear Compton scattering. An important point is that 
this theta function is generated from terms in both the 5 
function and principal value part of the photon propaga- 
tor and so the principal value part also plays a key role 
in the calculation of the two-step process involving a real 
photon. Even in the remaining one-step terms in Eq. ([7]) 

that comprise a cross-term and a direct-term it 
transpires that 



X«(r = 0)+^(r = 0)^0 



(12) 



So as the two-step term contains a contribution from the 
principal value part of the propagator so does the one- 
step term contain a non-zero on-shell contribution. Over- 
all causality can be seen to be observed without having 
to split up the propagator into 8 function and pricipal 
value par ts, by calculating the r-integral early on in Sfi 
using [23[ 



dr ■ 



1 



is 



3 *r(v.-vv) = -2m6(<p y -<p x ). (13) 



III. TWO-STEP PROCESS 

From the six original outgoing momentum integrals, 
due to the symmetry in the 1- (electric-field) direction, 
four integrals remain. As we neglect mixing between di- 
rect and exchange terms, each integral in the 2-direction 
can be factorised into a Compton-scattcring vertex part 
and a pair-creation vertex part. Some useful Airy in- 
tegrals that were derived from existing results in the 
literature are given in Sec. [B] which facilitated the 
p\ and p\ integrations, giving for P^ 2 > a final double- 
integral which can be written as a product of a spacctimc- 
dependent phase length squared and a dynamical part 
dependent on rclativistic and gauge-invariant quantities 
Xj = Xe{p° 3 - Pj), Xe = E/E CI and E cr = m 2 /e, 

= (ZAipflW/2 
j(2) = £^ f d X2dX3 0{xi -Xi- Xz) -4 (2) ( 14 ) 



(xi - X2Y 



where I (2) = I (2) (xi) and A {2) = -4 (2) (xi,X2,Xs) given 
by: 



A 2 



2/3 
M 2 



= Aii 
+a 3 Aii 



Aii 



2/3 
Ms 



■ a 2 Ai' 



2/3 

M 2 



Aii 



2/3 
M 3 



2/3" 
M 2 


Ai' 


2/3" 
M 3 


+ a 4 Ai' 


2/3" 
M 2 


Ai' 


2/3" 
M 3 



"2 

«:-! 



= M2 /3 (Xi+xMxi-X2) 
- = -M3 /3 [(Xi - X2 - X3) 2 + xl]/{Xi - X2) 
a 4 = -(M2Ms) 1/3 [xi - 2 X?(X2 + Xs) + X1X2MX 
-X2X3 + xl) + X?(2X2 + X2X3 + 2x1) + 
Xl(xi + 2X2X3 + 2x1)] /(X1-X2) 2 



M2 = 



Xl ~ X2 

X1X2 



M3 = 



Xi - X2 



(Xl - X2 - X3)X3' 



(15) 



(16) 



where Ai is the Airy function (24| , Ai' , its differential and 
Aii(x) = Ai(y)dy. The phase formation length can 
be written in a variety of ways: 



mxEAip _ L 
^0 ~ Tl 



(17) 



where L is the extent of the external field and = A/xe 
the formation length, where A = 1/m is the reduced 
Compton wavelength. It is quite logical that the for- 
mation length should be of this form when one considers 
over what extent the work done by the field is sufficient to 
create a pair L* = m/eE = A/xe- When an arbitrary ex- 
ternal field with phase ip and frequency x is taken as ap- 
proximately constant over the formation length one can 
apply the locally-constant-field approximation to P^ : 



P {2) (<P,<Po) 



dip x / d<p y \ dx2 dx3 



dl^ 2 



dX2dX: 



■(18) 



where a = (tuxe/* ) , Xi = Xi(<Po), X2 = X2 (¥>*), 
X3 = XziVy) an( i the phase integrals in Eq. (I18p arc 
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typical of other product approximations e.g. [25l. |26|| . 

Further support that the splitting off of the two step- 
process is unambiguous in a constant crossed field is 
found upon comparison with the "product approach" of 
integrating the tree-level processes of non-linear Comp- 
ton scattering (quantities denoted by subscript 7 ) and 
photon-seeded pair creation (subscript e ) over the inter- 
mediate photon lightfront momenta and summing over 
its polarisation, I. The differential rates for the sub- 
process are given by [22] : 




le-7 
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i 

XkZe 



Ai'(z 7 ) + Aii(z 7 



where z~, 



(Xfc/Xi(xi ~ Xk)) 



Ai'(z e ) + Aii(z e )j, 
(19) 

2/3 = ,r *- = 



(Xk/X3(Xk - X3)) 2/3 = /^3 /3 , and ± refers to transverse 
polarisations I = 1,2. The probability in the product 
approach is P ie = (^/S.tp) 2 X ie /2 where 



lA f xi dX y , 

-Lye = o / > / d Xk r, ' -Le,l- 

2 l~i Jo 9 *k 



(20) 



FIG. 2. The differential rate of the two-step process 
d 2 l (2) /dx'd X - for incoming fermion quantum non-linearity 
parameter a) xi — 1 an< l b) Xi = 10. As xi increases above 
1, the probability becomes skewed around large \ an d small 
X— 



The width of the distribution in p' 2 and p 2 _ can be esti- 
mated using the result that bremsstrahlung from an elec- 
tron is emitted in a cone of radius ~ 1/7 [28| . Assuming 
7 > 1, and that the incoming electron collides head-on 
with the external field wavevector, the magnitude of the 
transverse co-ordinate p^ ov cr the 3- co-ordinate of ini- 
tial electron momentum must be approximately equal to 
this angle, i.e. p' 2 jp\ « 1/7. Using the approximation 
X ~ 2"fXE, it then follows that p' 2 and hence p 2 _ must be 
approximately of the order of unity. 




By comparison of Eq. (TT9j) and Eq. (fToJI , one can observe 
that X^ = X~ e . This agrees with previous analyses 
performed for the total trident rate in a constant crossed 
field by analysing the polarisation operator correction 
to the Volkov propagator in a constant crossed field 

The differential rate d 2 I^ 2 ' /dx'dx- {x'> X- and X+ 
refer to the scattered electron, created electron and 
positron \ factors respectively) is plotted in Fig. [2] for 
incoming quantum non-linearity parameter xi = 1 and 
Xi = 10. It was found that for xi < lj the differential 
rate is symmetric in x' and x_, becoming more sharply 
peaked around the maximum at x' — X- = X+ — Xi/3 
as the initial x is shared equally among the products of 
the reaction. However when the highly quantum non- 
linear regime xi ^ 1 is entered, x' ~~ ► Xi and X± — ^ 1 
so the most probable scenario for a highly-relativistic 
fermion seed that creates a pair is that it remains 
highly-relativistic and can seed further generations of 
pairs in a field-free cascade (this point has recently been 
explored in [26|). 

By plotting the differential rate in the 2-momentum 
component of the outgoing electrons d 2 X^> / dp' 2 dp 2 _ in 
Fig. [3] it can be seen that at high xi, a beam of electrons 
colliding head-on with the external field wavevector is 
split into two connected phase-space regions, denoted by 
the two lobes in Fig. [3)3. Using transverse momentum 
conservation p\ —p l2 —p 2 _—p\ = 0, it can be shown that 
the distribution in electron and positron momenta is sim- where Xj — (Xj 1 -\-Xj2)/2 for j € {7, e}. The relative 
ilarly split, with the line of symmetry along p 2 , = — p_. difference between the total unpolarised and polarised 



FIG. 3. The scaled differential rate for the two-step pro- 
cess 10 s d 2 I^ 2 ' j dp' y dp y _ for incoming fermion quantum non- 
linearity parameter a) xi — 1 an d b) xi = 10- F° r higher 
Xi, the electrons are split into two ever-distincter momentum 
regions. 

The dependency of the total two-step probability on 
Xi is plotted in Fig. [5J in which 1^ is compared to the 
product approach using unpolarised tree-level differential 
rates 



■"ye 



dxk X e — y - 
°Xk 



(21) 
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rate was found to be around w 12% for 1 < xi < 10 3 - In 
addition, the dashed line in Fig. [4] plots an approximate 
formula derived by Baier, Katkov and Strakhovenko in 
Q, adapted to the current notation as: 



r(2) 



3a 2 



Xi 



^~T6^ l0g { 1 + W ° ~ (! + 0-56x1 +0.13 X f)s 

(22) 

which was found to agree with 1^ to within around 4% 
for 1 < xi < 10 3 . 
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FIG. 4. The two-step rate 2^ 2 ' is compared to an approxima- 

( 2) 

tion from the literature T^ ks and the product approximation 
of an integration over the lightfront photon momentum using 
unpolarised sub-processes I 7e . 



IV. ONE-STEP PROCESS 

P (i) 

was evaluated as a five-dimensional numerical in- 
tegral of the form 



(23) 



j(l) = [ d X2dX3dpldpldv6(xi - X2 - Xjj) 

TTXl 7 (Xl - X2) 2 

where 

B«(u) = |.4 (1) (t')--4 (1) (0)| 2 + 

2 Re (0) (a (1) * (v) ~A {1) * (0)) ] , (24) 

where = (u) = B {1) (v) + B {1) (-v) and A (1) (v) = 
•4. ( v i Xii X2, X3?P2i-Pi) are functions containing prod- 
ucts of Airy functions depending on the combination 

0!,3) 2 + "2,3 with: 



2/3 , 
M 2 3 ± 



1/3 
X1M23 



(25) 



where we note Vj(v = 0) = (fij/2) 2 ^ 3 , and the integral in 
v = 2xi>f°r/mx_E is between and 00. A check of the 



numerical expression was partly provided by the defini- 
tion 

h V \ dpj (0, xi , X2 , X3 , vl , vl ) = A {2} (Xi , X2 , X3 ) , 

(26) 

where A^ was given in Eq. (fTS")) . 

As noted in the derivation, the one-step term can be 
written as the sum of a purely one-step and a one-step- 
two-step interference term: 1^ = + \ where it 
was seen that X^. < 0. To investigate this point, we 
plot in Fig. [5] the differential rate d 2 I^/dx'dx-, in 
which one can clearly observe negative regions. As long 
as the total differential rate remains positive, this nega- 
tivity originates from the artificial splitting into one- and 
two-step terms. The one-step process in these momen- 
tum regions is not directly detectable but rather serves 
to suppress the overall rate. As xi increases, the sup- 
pression is reduced and the one-step probability becomes 
positive. 
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FIG. 5. A plot of the differential rate of the one-step process 
in x an d X- f° r a ) Xi = 1 an d b) Xi — 10. The one-step rate 
is entirely negative in a). 

To investigate the different dynamics of the one- 
step process, in Fig. [6] we plot the differential rate 
<9 2 zW jdp a dp 2 _^ in which one can clearly observe that 
the regions of maximum amplitude are negative. Upon 
comparison with the two-step differential rate in Fig. 
[3l we found that the purely one-step term 1^ did not 
split into two lobes as early with increasing xi an d the 
line of symmetry in the distribution was at a reduced 
angle. We consider whether the transverse momentum 
distribution can be used to measure the one-step process 
in the Sec. IVAl 

The total probability of the one-step process is plotted 
in Fig. [7] We note that the total one-step rate becomes 
positive for xi ^ 20. This is a consequence of the 
positive, purely one-step rate, which increases with xi 
becoming dominant over the negative cross term between 
one-step and two-step processes, in which the two-step 
rate decreases for xi > 20. It is somehow intuitive 
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FIG. 6. A plot of the differential rate of the one-step process 
in the 2-components of the outgoing electrons, multiplied by 
10 s when a) %i = 1 and b) xi = 10. 
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that the virtual process increases with xi- From the 
bandwidth theorem applied to the uncertainty relation, 
A£ AT ~ ft, the spread of virtual energies available 
A£ is larger when the interaction time AT is shorter, 
which here corresponds to the time taken to traverse the 
formation length T*. Since AT oc 1/ui, where u\ is the 
velocity of the incoming electron, which increases with 
Xi, the virtual photons can therefore more easily fulfil 
energy-momentum conservation at both vertices. The 
two-step process on the other hand, docs not benefit 
from this scaling. We discuss the consequences of this in 
Sec. EO 

From Fig. [7] we note that the maximum suppression of 
the two-step rate due to a negative one-step rate occurs 
at xi ~ 10. Since the cross-term is essentially the overlap 
between the one-step and two-step process, and since this 
grows with xij the suppression does not appear to be 
due to the increasing overlap of the two decay channels 
as in e.g. the LPM (Landau-Pomcranchuk-Migdal) effect 
(recently studied in the combination of laser and atomic 
fields [29|). What can be noticed, is that the probability 
for the one-step process becomes only positive when, over 
a formation length T», regardless of the extent of the 
field, the purely one-step process becomes more probable 
than the two-step process. 

The total rate for the one-step process can also be ver- 
ified in part by using existing asymptotic approximations 
in the literature, for example from Q one has: 

(i) -a 2 1 2xi , . 

^~^2-V^T e ' Xl<<1 (27) 

the error for which remains less than 20% for xi ^ 0.1, 
as plotted in Fig. [5] 

A test that was used to approximate the one-step 
background in E-144 [I, H[, is the Wcizsacker- Williams 
(WW) approximation P ww [H, [3(J , which is applicable 
to the purely virtual part, . This approximation sub- 
stitutes the virtual photon spectrum of a seed charged 



FIG. 7. Various parts of the one-step probability. Dotted 
lines represent suppression, and for xi ^ 20, the total one- 
step probability becomes positive (solid line). The purely di- 
rect term T^ 1 ^ is also compared with the Weizsacker- Williams 
approximation Xwi • 
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FIG. 8. The agreement between the one-step rate and the 
asymptotic expression T^kl fr° m 01 f° r Xi ^ 1- 

particle with an equivalent real photon spectrum (by 
neglecting non-transverse polarisation components) and 
then assumes the different real photon frequency com- 
ponents can be summed incoherently with the photon- 
seeded pair-creation probability. One can achieve a simi- 
lar result by performing the spin trace over just the trans- 
verse spacetime indices, followed by a limiting procedure 
k 2 — > in the resulting quantum amplitudes [311. This 
should be a good approximation to the trident process, 
the closer the intermediate virtual photon is to being real, 
i.e. the smaller k 2 becomes with respect to the electron 
mass. We note: 

k 2 = (P2 -Pi- tk) 2 (28) 
= {P2-Pif -2rx{p 2 ~pi), (29) 

therefore the closer the scattered electron momentum is 
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to the original one, i.e. the smaller the electron recoil, 
the smaller fc 2 is, independent of r. For the two-step case, 
the integration of the product approximation shows [26[ 
that in general \k <C Xi and X2 ~ Xi when xi ^ 1- It 
then follows that when xi 3> 1, the Weizsacker- Williams 
method can be used to approximate to the directly vir- 
tual part of the probability, P^ . One can modify the 
Weizsacker- Williams approximation for pair-creation via 
bremsstrahlung [22| to the present case 



Pe(Xk), (30) 



where C = j E + 1/2 - In 2 0.384, lE = 0.57721 ... is 
the Euler constant. The comparison is made in Fig. 
where we have taken P e = £AtpZ e and i^w = £A</?z4w- 
For values of xi ^ 2, the WW approximation becomes 
quickly worse than 10%, seeming to support the de- 
composition of the one-step process into purely-virtual 
and cross-term parts (although the error remains of the 
order of the polarisation error in the 2-step process, 
of around 10%). However, despite the accuracy when 
compared to the purely virtual one-step process, the 
WW approach does not take into account the cross-term 
between one- and two-step processes and therefore 
even fails to indicate the region of measurability of the 
one-step process or its suppressive effect on the total 
rate. One might conjecture that the WW rate could be 
useful when xi increases above 100, but we have not 
investigated this here as doubt has been cast on whether 
the perturbative expansion in final particles is valid at 
such large value of xi [Hf ■ 



critical at around xi ~ 10- One might argue that since 
L/ < 2, the 2^°) term must be taken into account and 
possibly cancels out this negativity. However, in Fig. [5Jd 
we plot the maximum ratio of L/L* required such that 
also the differential rate of the calculated terms in p' 2 and 
p 2 _ (solid line) and in x' and X- (dashed line) remains 
positive, which is seen to grow with xi into the L/L* > 2 
region {L/L* « 60 was found for xi = 10 3 , although it is 
unknown whether the current method is at all applicable 
at such high factors [HJ). Therefore the I^ a > term is 
unlikely to be fundamental to the discussion. 




FIG. 9. Plot a) is the sum of the total two- and one-step rates 
(multiplied by 10 6 ) as a function of incoming non- linear quan- 
tum parameter xi an d external field dimension L/L*. Plot 
b) is the maximum value of L/L* such that the differential 
rate in p' 2 and p 2 _ (solid line) and x' an d X- (dashed line) 
remain positive. 



V. TOTAL TRIDENT PROCESS RATE 

Comparison of the one- and two-step processes de- 
pends not only xi but also on the extent of the external 
field, and will allow for an analysis of how small the ex- 
ternal field can be taken without substantially violating 
the assumptions in the derivation. In previous sections, 
it was seen that for small xi> the total one-step process 
becomes negative. Since the total rate must remain posi- 
tive, this allows for a condition on the minimum allowable 
dimension of external field. To make the discussion more 
transparent let us rewrite the phase factors in terms of 
formation lengths recalling t;A.tp = L/L*, L» = X/xe, so 
that the total probability is of the form: 

P= I + (31) 

where 1^ corresponds to the neglected interference be- 
tween exchange terms. By regarding L/L* as a separate 
variable, we plot the dependency of the calculated terms 
in P on both xi and L/L* in Fig. [HI From the plot of the 
total rate in Fig. it seems that L/L* < 2 is required 
before the total rate becomes negative, and this is most 



The origin of this negativity is most certainly the as- 
sumptions made in deriving the constant crossed field 
rates. Most notably, the infinite integral over the 
external-field phase performed to generate the charac- 
teristic Airy functions must be modified to include the 
finite duration. These results imply a useful additional 
constraint on the validity of the constant crossed field 
approximation, namely 

£Atp > 1. (32) 

This is independent of the condition £ 3> 1, which is re- 
quired to justify the limit on the external- field frequency 
x {) — > in the expression for a general background 
because £A<p is independent of x°. 



A. Measurability 

Here we are discussing the measurability of electron- 
seeded pair creation in the collision of a laser and elec- 
tron beam in the non-perturbative and highly non-linear 
regime (£ 3> 1) in contrast to the weakly nonlinear regime 
(£ ~ 1) discussed in (BJ in a curtailed plane- wave back- 
ground. The two-step process is measurable by tak- 
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FIG. 10. In plot a) for xi — 10 is the ratio of the one-step 
to the two-step process, where in the empty regions, the ratio 
becomes much larger than the maximum on the colour scale. 
In plot b), the percentage r(xi) of the measurable one-step 
signal to the total signal in the p 2 _-p\ plane, with increasing 
I//L*, is plotted. 



ing the similar parameters to the SLAC E-144 experi- 
ment, but increasing the intensity of the laser used from 
1.3 x 10 18 Wcm~ 2 to 10 22 Wcm -2 and reducing the pulse 
duration from r = 1.6 ps to r = 10 fs allowing the en- 
ergy and number of initial electrons to be reduced to 10 9 , 
3.5 GeV electrons. Assuming the slight angle makes only 
a minor difference to the two-step probability, (xi) = 1, 
£ = 30, LjL % = 1100, where (•} refers to an average over 
oscillations of the external field, the initial beam of elec- 
trons would create of the order of 10 6 pairs. We stress 
that Compton scattering and coherence effects have been 
neglected. The one-step process in this regime seems 
much more difficult to separate in an experiment. One 
possibility is to measure the transverse momenta of final 
state positrons (the electron momentum distribution will 
likely be washed out by double-photon Compton scat- 
tering of the initial electron). Due to the slightly wider 
momentum spectrum in the 2-direction for positrons cre- 
ated via the one-step process compared to the two-step 
process, plotted in Fig. 110a . with a judicious momentum 
"cut", the effect of the one-step process could be mea- 
sured. In Fig. [TUb the ratio r(xi) of one-step to two-step 
created positrons outside of the area < m, in regions 
of the detector where the former is at least 50% the latter, 
is plotted as a function of the field dimension L/L*. The 
one-step process is then most measurable for a head-on 
collision when the ^-factor is increased but xe held rela- 
tively low. A beam of 10 9 , 500 GeV electrons combined 
with a single-cycle 10 21 Wcm -2 laser pulse ((xi) = 20, 
£; = 10, L/L* « 60) would produce again of the or- 
der of 10 4 positrons via the one-step process in these 
measurable regions from initial seed electrons. Although 
further Compton-scattcring of the positrons could take 
place, this three- rather than a two- vertex dressed pro- 
cess is less probable. 



VI. DISCUSSION 

A fundamental difference between the two- and one- 
step processes is that in the two-step process, the real 
photon can propagate for an arbitrarily long phase 
length before decaying into a pair, whereas in the 
one-step process, the photon's range is limited by the 
uncertainty relation. So when the sub-processes are 
not very probable, for example at low Xi one would 
expect the two-step process to dominate. If the range 
of the real photon were curtailed to a length of the 
order of the formation length, the length associated with 
the one-step process, this advantage of the two-step 
process would be lost and the two processes should be 
comparable, as we indeed find. On the other hand, 
considering the uncertainty relation in energy and time, 
the shorter the interaction duration, the larger the 
bandwidth of frequencies available to fulfill energy- 
momentum conservation at the vertices of one-step pair 
creation. Therefore, for electrons incident with a high 
X parameter, the probability of the one-step process is 
enhanced compared to the two-step one, leading to its 
dominance for external field dimensions of the order of 
the formation length and large x- 

An important point is the overlap of the two pro- 
duction channels. Assuming that the constant crossed 
field is still valid (£A<y9 ;§> 1), instead of considering 
the suppression of the two-step through the one-step 
process, one could instead consider this phenomenon 
the other way round. The one-step process is es- 
sentially suppressed by the more probable two-step 
process at low Xi hut when x is large enough that 
the one-step channel becomes more probable over its 
formation length, the one-step process probability 
becomes positive and emerges as a separate channel 
on its own. Moreover, this threshold is independent of 
the extent of the external field and so this phenomenon 
is likely present in more complicated backgrounds as well. 

Although we have only calculated a two-step process, 
the question presents itself whether the constant crossed 
field approximation can be used to model long chains of 
processes as is typically used in computer codes. In this 
case for n vertices, assuming the Poisson-like dependence 
on the external field's spatial extent (and duration) con- 
tinues to higher orders, one could have a probability of 
the form: 



P 



1 



L 



XW+... + ^-lW+tf 0) , (33) 



where represents a process with j — 1 on-shell 

propagators. Assuming the dynamical factors 1^' 
are of the same order, the potentially largest negative 
contribution arises from the xOO-ZO" 1 ) cross-term, 
which acquires an extra factor n. Therefore one could 
surmise that the constant-crossed-field approximation is 
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also valid for an n-vertex process as long as £Aip ^> n. 

The two-step process was measured in the weakly 
non-linear (multi-photon) regime (£ = 0.3) in the E- 
144 experiment. We calculated the approximate num- 
ber of pairs created in the highly non-linear and non- 
perturbative regime (£ ^> 1, % > 1). It was found 
that if the laser intensity could be updated to higher, 
currently available intensities (2 x 10 22 Wcm -2 was al- 
ready achieved in 2008 [20]), the particle beam can be 
allowed to be even less relativistic and the two-step pro- 
cess measurable. More difficult is separating the one-step 
mechanism for which it is crucial that the field extent is 
not much larger than the formation length otherwise the 
background from the two-step channel is too large. How- 
ever, to utilise the slightly different transverse positron 
momentum distribution, the x factor must remain high 
X ^ 10. Here the same order of number of positrons was 
predicted to be produced in "measurable" regions for a 
single-cycle 10 21 Wcm -2 pulse but a 500 GeV electron 
beam. 



VII. CONCLUSION 

The trident process in a constant crossed field must be 
considered in its entirety, being separable into two- and 
one- step processes rather than real and virtual parts, 
which were both seen to contribute to the one-step pro- 
cess. The two-step process was found to agree exactly 
with an integration of the average of polarised tree-level 
processes over hghtfront momenta. The one-step pro- 
cess was found to be effectively suppressed for xi ^5 20 
due to the larger probability of the two-step process over 
the formation length. In the highly non-linear and non- 
perturbative regime 3> 1, x > 1), the two-step process 
was shown to be measurable for electron beams even less 
energetic than in the E-144 experiment, as long as the 
intensity of the laser is updated to around 10 22 Wcm -2 . 
For field dimensions not orders of magnitude larger than 
the formation length, it was shown that the one-step pro- 
cess could in principle be separated from the two-step 
process using the wider positron transverse momentum 
spectrum when a single-cycle 10 21 Wcm -2 pulse collides 
with a 500 GeV electron beam 
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Appendix A: Further detail on derivation 
1. Definitions 

Here we define objects used in the manuscript and fur- 
ther calculation. The Volkov states are Hal: 



ip r (p) 
S(p) 



l + eph u r (p) piS(p) 
2xpi y / 2p°V 



y f 2 P W 
1 _ e ^ 



2 Hp 



2xp . 

v r {p) 



y/2pW 



iS(p) 



iS(-p) 



-px 



i(pA\g) _ e 2 A 2 [<f>] 
2(*p) 



Hp 



(Al) 



where = 7 M >f M , 7^ are the gamma-matrices, 
u r (v r ) are free-electron (-positron) spinors sat- 
isfying J2l=i u r p(p)u rr7 (p) = (f + m) pa /2m, 

J2l=i v r P (p)vra{p) = {j/> - m) pcr /2m, u = u^j° and 
the remaining symbols are as described in the paper. 
The photon propagator is: 



d 4 k 4ng 



ik(x-y) 



(2tt) 4 k 2 + is 
The expansion of the vertices in Fourier modes is: 



(A2) 



fly 



A" (si 



(A3) 

^e !r ^ 2 M 7 ^iM (A4) 
(A5) 

dip e^^3M7^4 + M(A6) 



where we have used the shorthand ipi = ip(pi) with spinor 
indices suppressed and ipj(<p) are the Volkov states with 
Fourier terms e PjX removed. 



2. Derivation of rate expression 

Beginning from the expression for the scattering ma- 
trix: 



S/i = a I d 4 xd 4 y ip 2 (x)'y> i ip 1 (x)G llv (x - y)tp 3 (y)j v ipi (y) 
-(ft^R), (A7) 
= Sft-%7i- (A8) 

Using the definitions in Eqs. (|AHIA6j) . we can rewrite 
Eq. (fMI as: 

^ = ^[d 4 xd 4 y ^-drds e^^ ElMl , 
1 ttJ y (2tt) 4 k 2 +ie 

(A9) 
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where Ilr = k — 5pr, Spr = Pi + rx — P2 and I1a = 
— k — Spa, Spa = sx — P3 — Pi- Integration of Eq. (|A9|) 
over x and y gives: 



Sfi = a Jd 4 kdrds S(k - Sp r ) 



xS(k + $pa)- 



k 2 + is 



and over k gives: 



Sfi - 



(27r) 4 a 



drds S^(Ap - (r + s)x) 
r^(r)A M (s) 



fc' 2 



fc'=fe' 



(A10) 



(All) 



where k^ — Sp + rx, dp 
Pi — p\. We notice: 

1 



k' 2 + 



is 



Pi - P2 and Ap = p 2 + p 3 
1 pttSp)- 1 



k'=s P +r>t (Sp) 2 + 2rxSp + is 



r — r* + is 
(A12) 



where we have defined = — (5p) 2 / (2x5p). With a co- 
ordinate transformation r — > r + r* we have: 



(27r) 3 a fdrds 
>r<5p J ?* + is 



<5 (4) (Ap- (r + r* + s)x) 
r"(r + r.)A M (s), (A13) 



In order to evaluate the delta functions, we switch at this 
point to lightcone co-ordinates. For spatial co-ordinates 
we define x = (x° ± x 3 ), x± = (#1,2:2), whereas for 
momenta, p ± = (p°±p 3 )/2, p± = (p 1 ,p 2 ). We also define 
a co-ordinate system and specify a constant crossed field 
h = x°(l, 0, 0, 1) A^((p) = a^ip, = (E/x°)(0, 1, 0, 0), 
xa = >c 2 = 0, so that xx = x°(x° — x 3 ) = x + x~ . 

In forming the probability, we must square the scatter- 
ing matrix. Let us concentrate on |S/i| 2 as the steps for 
other contributions are similar. When Eq. (| Al 1|) is mod- 
squared, one has, for some function / = f(r,s,r',s') £ 



\s\ 2 
\s\ 2 
\s\ 2 
\s\ 2 
\s\ 2 



dr dr' ds ds' f <5 (4) [Ap - (r + s)x] S (4) [Ap - (/ + s')x] 

dr dr' ds ds' f S {4) [Ap - (r + s)x\ S {4) [(r + s - (r + s'))x] = 

drdr'dsds' fS^[Ap-(r + s)x] 8^ [(r + s- (r' + s'))x] 6 Jf + S - t~, + 8 '}1 

o(r + s — (r + s'j) 



/ dr dr' ds ds' f 5^ [Ap - (r + s)x] , 
J (2tt) 6 A<p + 



5(r + s- (r' + s')) 



VT 



(x ) 4 (2tt) 3 A^ + 



<5 (2) (Ap- L )(5(Ap-) J drdr' f(s = Ap-r,s' = Ap-r'), 



(A14) 
(A15) 
(A16) 
(A17) 
(A18) 



r 



where we have defined a formation phase length Aip + 



S(r + s — (r + s )) 



r-^s—r'-^s' 2,71 



where 



lx=0 








/ 2tt 


x=0 



(A19) 



(A20) 



At this point, since we wish to form probabilities and 
not rates, we invoke the relation Tjp\ = Aip + /xpi [Ioj |. 
so that, combining the arguments leading to Eq. (|A13|) 



and Eq. (|Al8)l . we then have: 
-> 2 (27r) 3 a 2 Vp\ 



Z(->,->) 



(x6p) 2 (x°) 4 (x Pl ) 



S^iAp^SiAp-)!^,- 



dr 



rn(r + r*)A^(s*-r) 



r + is 



(A21) 



where we have defined s* = Ap + — r*, which can be 
shown to be equal to: 



(j>2 +P3 +Pi) 2 - m 2 
2pix/x° 



(A22) 



We note that in order to evaluate the light-cone co- 
ordinate delta functions occurring in Eq. (|A21[) from 
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a Cartesian integral, one can use: 
d 3 p 



2 / f / ' (/M ; J 2p 



d 2 pidp 



where (?(•) is the Hcaviside step function. 

The probability ~P, using the expression P 

(1/2) nUlvf^Pi/i 2 *) 3 }* 1 l^l 2 ' is thcn s iven h r- 



n 



d 2 pf dpj 0(p7)trX(-»,-») 



(2x0)«(xp 1 ) .£J- 3 y (2tt)3 p - p7(p- 



(A24) 



where the instruction nn means that all normalisations 
of the form (2Vp^)~ 1 / 2 for j € {1,2,3,4} have been re- 
moved and the integral in d 3 p^ has already been per- 
formed (to account for the degeneracy of outgoing states 
the total probability P requires an extra factor 1/2 as 
explained in the main text). 



3. Vertex functions 

We can rewrite the vertex functions Eqs. (|A4|) and 
(|A6|) in a way that allows them to be easily evaluated by 
separating integrals from trace products. Concentrating 
first on r M (r): 



1*0) = dip 



u a2 (gg) 



s4(<p)# 



i 



2xrpi 

where we have introduced 



2xp 2 

u <?l (Pi) t(ry+c 2 y 2 +c 3 y 3 ) 



A25) 



efp 2 a pia\ 

c 2 = o ( ! c 3 

2 V Xp2 Xpi / 



e 2 a 2 / 1 1\ 

6 V xp 2 xpi ' 

(A26) 



Now as A M = a M (p, we can rewrite Eq. (|A25jl as 

U<t 2 (P2) 



r*(r) 



n/^P 7 



4xp2>«Pi 



(A27) 



C„(r,c 2 ,c 3 ) = / ^"- 1 e i ^ +c ^ +C3 ^ >.(A28) 

J — oo 

By shifting the </? co-ordinate ip ip — C2/3C3, one can 



show: 

Ci = 6Ai(/i 2 / 3 ) 



c 2 = - 

C 3 = b 



(3^ Ai ' < " 2,a) + S A1( '' 2/a) 



(A29) 
(A30) 



2/3 



a Ai <" 2/3 ) 



2ic 2 



Ai'(^ 2 / 3 ) 



(A31) 



2- 



(3c 3 )V3 



V = 



rc 2 2c 3 , 



3c 3 27c 2 .' 



,2/3 _ 



c|/3c 3 



(3c 3 )V3 



(A32) 



where Ai is the Airy-function defined in Eq. (|B12[) with 
normalisation N = tt. 

The calculation of A AI (s) proceeds in a similar way. 
From Eq. (|A6|) we have: 



A M («) 



dip J " g3 ^ [1 + e ^)^ 



2^ P3 y 



7p 



1 _ e^4((^) i ^ 4 (p4) pi(ay+ 4 y 3 +4y 3 ) 



2^P4 J V2p°V 
where we have introduced: 



V33) 



, e/p 3 a p 4 a\ , 

c 2 = o ' c 3 

2 V >ip 3 >rp4 / 



e 2 a 2 / 1 



6 V xp 3 xpi ) 

(A34) 



Then: 



M*) 



?V 3 (P3 



-£>3 



2 \ xp3 xp4 / 



4>fp 3 >£p4 



«g 4 (P4) 



(A35) 



D„ = C n (r -> s, c 2 ->■ c 2 , c 3 ->■ c 3 ). (A36) 



4. Fermion trace 

The trace to be evaluated comprises the trace of each 
exchange term mod-squared plus interference terms. If 
we define the objects: 



AI, = 7°M+7 ; 



2m 



(A37) 
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then for each exchange term squared, the trace is of the 
form: 



►) = ^ti\u a2 C^(p2,Pi,r)ii ai u (y3 D t _ l (p 3 ,p 4 ,r)v a 



u a2 C' y (p2,Pi,r')u ai u a3 D„(p 3 ,p 4 ,r')v 



-tr 



£+C^{ P2 ,p u r)£+C^(p 2 ,p 1 ,r f ) 



.tr 



£ t D fj.(P3,P4,r)£ 4 Dl(p 3 ,p 4 ,r') 



(A38) 



where from the definition of I(— >, — >) Eq. (|A21|1 . C M and 
are factors of Cj and Dj, j 6 {1, 2, 3}, multiplied by 
the combinations of gamma matrices occurring in Eqs. 
(IA32I) and (1A35|) . integrated over r and r' variables. Fol- 
lowing similar steps, for each interference term it is of the 
form: 



) =zl tr u <ys C>l {P3;Pi, , r)u <7l u a2 D ll ,(p2,Pi,r)v a 



u rJ2 C / (p2,Pi,r')u ai u a3 D l ,(p 3 ,p4,r / )vo 



-tr 



£+C^(p2,p 1 ,r')£ 2 tD^p2, Pi ,r) 



£ 4 Dt(p 3 ,p 4 ,r')£+C^p 3 , Pl/ , 



(A39) 



These traces were performed with the package Feyncalc 



5. Complex phase factor 

It can be seen from the definitions of the Airy integrals 
resulting from the vertex factors, that an overall phase 
factor r](-, •) occurs in the traces •) (from the r\ factors 
in Eq. (TAl2"|) occurring in Eqs. §KM) and (|Al9]) ). If one 
squares the r integral, labelling the new co-ordinate r', 
they are of the form: 



v(->,->) = - 



C 2 (P2,Pl) _ c' 2 {p 3 ,P4) \ , _ n 

3c 3 (p 2 ,pi) 3c' 3 {p 3 ,p 4 )J 
(r + r„)c 2 {p2,Pi) . (r - s*)c 2 (p 3 ,p 4 ) 



3c 3 (p 2 ,Pi) 3c' 3 (p 3 ,p 4 ) 
(r' + r <f )c 2 {p 3n pi) _ (r' - s*)c 2 {p 2 ,Pi) 
3c 3 (p 3 ,pi) 3c' 3 (p 2 ,p 4 ) 

2 / cl(p 2 ,Pl) 4(P3,Pl) . c' 2 3 (p 3 ,p 4 ) 



27 \cl(p 2 ,pi) c§(p 3 ,pi) c' 3 2 {p 3l p 4 ) 

C' 2 3 (P2,P4) 
C' 3 (P2,P4) 

(A40) 

After some simplification, it can be seen that the non- 
interference terms have a relatively simple structure 



(where xe = E/E cr , E CI = m 2 /e): 
-x°(r-r') 



m 2 XE(p 1 ~P 2 ) 
+P{(P2 ~Ps )] 



— [pI(P3 -Pi ) +p\{Pi ~Pi ) 

(A41) 



6. Formation length of two-step and sub-processes 

If we imagine that the process takes place in a co- 
herence interval of finite duration, for a given incoming 
electron momentum p\ , these integrals will also be finite 
(otherwise the particles would have to be accelerated in- 
finitely quickly (see also [3)). We apply the following 
reasoning, which is standard for lower-order constant- 
crossed- field processes (see e.g. H3])- The phase of the 
modified Airy functions that occur at each vertex Eq. 
(|A28[) have a stationary point at: 



1 ± 



3rc 3 



P 



c 2 
3c 3 



(A42) 



Let us write ip* = —p(l ± Sip*), where p is the phase 
at which the process takes place (the saddle-point) and 
Sip* is the width. It is assumed that Sip* = because 
when this is not the case, the resulting Airy functions 
have no dependency on dynamical variables. Therefore, 
integration over p is equivalent to integration over the 
relevant part of the phase. 

The complex phase factor for the purely direct (and 
analogously for the purely exchange) term in a single r- 
intcgral is of the form: 



T)(- 



(A43) 



= (<Px ~ Vl)r, 

where ip* z is the saddle-point at co-ordinate z and for the 
purely direct term 



x° p\p x 



P\P 2 



m 2 XE 



<Pv = 



Pi ~Pi 

x° p\p 3 + p\ (Pi -P 2 )- PiPz 
m 2 XE p 4 - P2 



(A44) 



(A45) 



7. Isolation of the two-step process 

Concentrating on the non-exchange term — >) (an 
analogous calculation follows for !(•<—,•<—)), the integral 
over r is of the form: 



1 



J =^2 \ d P\ d P\ 



dr 



e i[fl{pl)-V*y(phpl)} r pU 



(r + is) 



(A46) 



with F(r) 6 C. Performing an integral substitution ip*± 
ip*±ip*, one can rewrite this as: 



J 



2Jir 2 



dip*, dip*_ 



dr 



F(r) 



(r + ie) 



, (A47) 
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where J = \d(<p+, Lp*_) / d(p\ : p\)\ ^ s the inverse Jacobian. 
In order to remain consistent, before integrating in the 
variable (p * , wc will first perform the principal value cal- 
culation. The order of integration is important as prin- 



ciple value an d Lp* + _ integrals do not commute (for ex- 
ample in Eq. (|A51j) . integration in a does not commute 
with the operation V). 



J 



1 



2Jtt 2 
1 

2Jn 
1 
J 
1 
n 



dtp* + d(p*_ 
! dip* + dip*_ 



-iirF(0)+P dr 



e lv ~ r F{r) 



-2mF(G)6(-ip*_) + dr e lv - 



_ r F(r)-F(0) 



2|F(0)| 2 / dipl dip*_ 6(-tp*_) + - dip 



F(0) / d V 



dr 



dr 



\F(r)-F(0)\< 



F*(r) + F*(-r) - 2F*(0) 



(A48) 
(A49) 

(A50) 



where c. c. stands for complex conjugate and in Eq. 
(|A49|) and Eq. (|A50|) respectively, we have used the re- 
sults (H: 



V 



dr 



in sgn(a) 



1 



(A51) 



d<p0{ip)e lipr = iV-+Tv6{r). (A52) 
r 



The first term in Eq. (|A50|) can be identified as the two- 
step process due to the two phase integrals that occur, 
justified in the main text. In addition, a Heaviside theta- 
function in — ip*_ occurs, which is a sign that causality is 
preserved insofar as pair-creation can only occur after 
Compton scattering in this setting. 



8. Justification for neglecting interference terms 
between direct and exchange parts 

For the non-interference terms, it has been shown in 
Eqs. (|A46HA50f how the simple nature of the exponential 
occurring in the vertex functions leads to a dependency 
on the external field phase in Aip. For the interference 
terms the integral over p\ 3 is of the form 



J (-+,<-)= \dp\dp\drdr' 



M^,^-)F(r)F*{r') 
(r + is)(r' — is) 



<A53) 



where the phase rj(— y,<— ) contains terms of the order 
P2,3: (P2.3) 3 ; pUpI) 2 and (P2) 2 Pl- Instead of generating 
an arbitrarily large (divergent in the strict sense) fac- 
tor of formation length Aip as for each exchange term 
squared, for this interference term, an Airy function in 
remaining particle momenta and r is generated, which, 
having positive or negative argument, will tend to reduce 
the value of the integral. As noted in the conclusions of 
the main text, if one demands that £Aip 3> 1 to prevent 



I 

regions of negative total probability from arising, then 
these interference terms can be safely neglected. 



9. Numerical evaluation of one-step integral 

Despite the five-dimensional integration of 2A 1 ) Eq. 
(|A54[) not being oscillatory, it is challenging to numer- 
ically evaluate 



o? f dx2dxzdpldpldv6(xi ~ X2 - Xs) # ( / 



(Xi ~ X2) 2 



^A54) 



As the evaluation of B^ 1 ' is computationally expensive, it 
is important to know the relevant bounds of the variables. 
In the main text, it was justified that p\ 3 f=a 1. We 
then take £>!>P§>£ [~ 4, 4]. From the arguments of the 
Airy function given in Eq. (|25|) . one would expect the 
maximum of in v to be of the order v « XiP-2,3- Let 
a = X2/X1 an d b — X3/X1 so that a, 6 £ [0, 1], then 



X1P2 = - - 1; 
a 



1 



XiM3 = 



1 



\ (A55) 



However, from studies of the approximated two-step 
process in [2(|, it seems that as Xi increases above 1, 
(1 — a) /a ~ 1/xi- Likewise, X3 was observed to remain 
approximately constant so that b — > l/xi> leading to 
XiA < 2,3 ~ Xi- Therefore, v e [0, lOxi] was chosen for the 
v integration, with the tail [lOxi, 00] evaluated in w with 
the conformal transformation w = tan -1 v. Although the 
function B is quite smooth in the X2 _ X3 plane, the largest 
contribution to the total integral originates from an ever- 
smaller region around a = 1, b = 0, with increasing xij 
making these points particularly costly to evaluate. To 
escape the triangular X2~X3 plane as given in Fig. 2 in 
the main text, one can substitute integration variables 
X2 -> X2/(xi - X3) and X3 -» Xa/Xi to achieve a square 
integration region between and 1. One can then easier 
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observe where the maxima lie in the integrand and eval- 
uate grids of points incorporating these. The resulting 
surface can then be interpolated and numerically inte- 
grated. For x > 1, The global maximum was found in an 




n 



6.4e-06 



4.0e-06 



■ 2.0e-06 

I 



0.93 0.94 0.95 0.96 0.97 0.98 0.99 1.00 

X2AX1 ~ Xs) 



FIG. 11. The differential of the dynamical part of the one-step 
rate dl w /dx2d%3 for X i = WO. 

ever smaller region X3/X1 ~> 0, X2/(Xi — X3) — > 1- For 
example, for the case X i = 100, as plotted in Fig. [Ill 
the maximum was centred around Xz/iXi ~ Xa) = 0.995, 
Xs/Xi = 0.005. The tests of accurate integration, other 
than variation of number of points and integration re- 
gion, were provided by comparison with asymptotic ex- 
pressions from the literature for the total rate and the 
Weizsacker- Williams approximation. 



Appendix B: Integrals of Airy functions 



where c is an arbitrary constant (integrals involving coef- 
ficients with odd powers of t are zero due to the functions 
being odd). With Iq, J and K being given in e.g. [3o| . 
from partial integration and the use of some primitives 
given in [361 ] . the following analytical results have also 
been verified numerically: 



h 
h 

h 
Jo 

K 2 



7T 

16JV 

7T 



Ai'(u) + cAii(u) 



4Nk 

7T 

~8Nk 
3tt 



^Ai(w) + -Ai'Ov) + c 2 Aii(z;) 

. 4 K 

Ai(v) 



16N 
n 



Ah(v) 
1 

K 



4Nk 

7T 

16N 



Ai'(u) + cAii(w) 

[3Ai'(v) +ckM 1 (v)] 
5 c 

—kM(v) Ai'(v) - c 2 Aii(w) 

4 K 



(B4) 

(B5) 

(B6) 
(B7) 
(B8) 

(B9) 
(B10) 

(511) 



v = kc, k = 2 2 / 3 , N is the normalisation factor occuring 
in the definition of the Airy function: 



1 f°° 

Ai(x) = — / dt cos (t 3 + xt) (B12) 



We give here a selection of Airy integrals that are use- 
ful in the derivation and are in part derived from other 
results in the literature. Let us define: 



Jin 
K 2n 



dt t 2n Ai 2 {t 2 +c) (Bl) 
dtt 2n Ai(t 2 + c)Ai'(t 2 +c) (B2) 
dt t 2n Ai l2 {t 2 + c), (B3) 



where Ai (x) = dAi(x)/dx and Aii is defined as: 



Aii(cc) = dt Ai(t + x). 
Jo 



(B13) 
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